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n this paper, we discuss improvements of the 
Suto et al. (2000) model, in the hght of recent the- 
oretical developments (new theoretical mass func- 
tions, a more accurate mass-temperature relation 
and an improved bias model) to predict the clus- 
tering properties of galaxy clusters and to obtain 
constraints on cosmological parameters. We re- 
derive the two-point correlation function of clus- 
ters of galaxies for OCDM and ACDM cosmolog- 
ical models, and we compare these results with 
the observed spatial correlation function for clus- 
ters in RASSl (ROSAT All-Sky Survey 1), and 
in XBACs (X-RAY Brighest Abeh-Type) samples. 
The comparison shows that the best agreement is 
obtained for the ACDM model with fim = 0.3. 
The values of the correlation length obtained, 
(ro ~ 28.2 ± 5.2h-i Mpc for ACDM), are larger 
than those found in the literature and compara- 
ble with the results found in Borgani, Plionis & 
Kolokotronis (1999). In order to study the possi- 
ble dependence of the clustering properties of the 
X-ray clusters on the observational characteristics 
defining the survey, we calculated the values of the 
correlation length rg in the catalogues where we 
vary the limiting X-ray flux S'lim- The result shows 
an increase of ro with Liim, and correlation lengths 
that are larger than in previous papers in litera- 
ture (e.g Moscardini, Matarrese & Mo 2001 (here- 
after MMM); Suto et al. 2000). These differences 
are due essentially to the different M-T, mass func- 
tion and bias model used in this paper. Then, we 
perform a maximum-likelihood analysis by com- 
paring the theoretical predictions to a set of ob- 
servational data in the X-ray band (RASSl Bright 
Sample, BCS (Rosat Brightest Cluster Sample), 
XBACs, REFLEX (ROSAT-ESO Flux Limited X- 
Ray Sample)), similarly to MMM. In the frame- 



work of cold dark matter models, we compute the 
constraints on cosmological parameters, such as 
the matter density fl^i, the contribution to density 
due to the cosmological constant, VIa, the power- 
spectrum shape parameter F and normalization 
cTg- If we fix r and as, at the values suggested 
by different observational datasets, we obtain (for 
flat cosmological models with varying cosmologi- 
cal constant Qqa = 1 — i^om) constraints on the 
matter density parameter: 0.25 < riom < 0.45 and 
0.23 < r^om < 0.52 at the 95.4 and 99.73 per 
cent levels, respectively, which is 20-30 % larger 
than the values obtained MMM. Leaving F, and 
ftmo, free for the flat model, the constraints for 
F are 0.1 < F < 0.14, while for the open model 
0.09 < F < 0.13. These values are smaller than 
those of MMM by about 20 - 30%. If we keep the 
values of J^a fixed, we obtain the constraints in the 
T — as plane. For the open model with f^om = 0.3 
the 2a region for F is 0.11-0.2 for as it is 0.7 and 
1.55. For the flat model with fiom = 0.3 the 2a 
region has 0.13 < F < 0.2 and 0.8 < as < 1.6 
The values of as obtained are larger than those of 
MMM by ~ 20%. If we allow the shape parame- 
ter to vary, we find that the clustering properties 
of clusters are almost independent of the matter 
density parameter and of the presence of a cosmo- 
logical constant, while they appear to be strongly 
dependent on the shape parameter. 

Key words, cosmology: theory - large scale struc- 
ture of universe - galaxies: formation 



1. Introduction 

X-ray studies of clusters of galaxies have provided 
a large amount of quantitative data for the study 
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of cosmology. The mass of a rich cluster is ap- 
proximately lO^^h~^M0. ^ This mass lies within 
a region of diameter ~ 20h~^Mpc and conse- 
quently the observations of clusters can provide in- 
formation on the mass distribution of the Universe 
on these scales. Furthermore, since rich clusters 
are rare objects, their properties are expected to 
be sensitive to the underlying mass density field 
from which they originated. Therefore, clusters of 
galaxies appear to be ideal tools for testing the- 
ories of structure formation as well as studying 
large-scale structure. 

X-ray catalogues of X-ray selected clusters 
arc now available from ROSAT: RASSl (De 
Grandi et al. 1999) (ROSAT All-Sky survey 1), 
BCS (EbcHng et al. 1998) (ROSAT Brightest 
Cluster Sample), XBACs (Ebehng et al. 1996) 
(X-ray brightest Abell Cluster Sample), REFLEX 
(Bohringer et al. 1998) (ROSAT-ESO Flux- 
Limited X-ray sample) and the volume covered 
by the samples is expected to increase through 
the X-ray satellites such asAstro-E, Chandra, and 
XMM. These data together with optical data have 
been used to compute the cluster number counts 
and the X-ray luminosity function, which have rel- 
evant cosmological implications. In particular, the 
analysis of the cluster abundance (also as a func- 
tion of rcdshift) has been used widely to provide 
estimates of the mass fluctuation amplitude and 
of the matter density parameter Q^, with sev- 
eral, often discrepant results (Kitayama & Siito 
1997: Mathiesen & Evrard 1998; Sadat, Blancliard 
& Oukbir 1998; Reichart et al. 1999a,b; Viana & 
Liddle 1999; Blanchard, Bartlctt & Sadat 1998; 
Eke et al. 1998; Bahcall, Fan & Con 1997; Fan, 
Bahcah & Cen 1997; Bahcall & Fan 1998; Donahue 
& Voit 1999; Borgani et al. 2001). An alternative 
approach to the abundance of clusters is based on 
the study of the spatial distribution of selected 
clusters. The standard statistical tools used with 
this aim are the (spatial and angular) two-point 
correlation function and the power-spectrum. 

The two-point correlation function is a funda- 
mental statistical test for the study of the clus- 
ter distribution and is relatively straightforward 
to determine from observational data. The spa- 
tial correlation function of galaxy clusters pro- 
vides an important cosmological test, as both the 
amplitude of the correlation function and its de- 
pendence upon the mean intercluster separation 
are determined by the underlying cosmological 
model. Like for the cluster abundance, discrepant 
results have been found for, e.g., the correlation 

^ h = Ho/(100kms"^Mpc-^), Ho being the Hubble 
constant at the current epoch (in the paper we adopt 
h = 0.65) (see Spergel et al. (2003), Tegmark et al. 
(2004)). 



length (Hauser & Peebles 1973; Bahcall & Soneira 
1983; Klypin & Kopylov 1983; BahcaU & Cen 
1992, Bahcall & West 1992; Efstathiou et al. 1992, 
Governato et al. 1999). 

As shown in some papers (Eke et al. 1998; 
Reichart et al. 1999a,b; Donahue & Voit 1999; 
Borgani et al. 2001; Del Popolo 2003), the rea- 
sons leading to the quoted discrepancies are not 
only connected to the observational data used, but 
other unknown systematic effects may be plagu- 
ing a large part of the quoted results (Reichart 
et al. 1999a,b; Eke et al. 1998; Donahue & Voit 
1999; Borgani et al. 2001). Systematic effects en- 
tering the quoted analyses are: 1) The inadequate 
approximation given by the Prcss-Schechter rela- 
tion (e.g., Bryan & Norman 1998). 2) Inadequacy 
in the structure formation as described by the 
spherical model leading to changes in the thresh- 
old parameter Sc (e.g., Governato et al. 1999). 
3) Inadequacy in the M-T relation obtained from 
the virial theorem (see Voit & Donahue 1998; 
Del Popolo 2002a). 4) Effects of cooling flows. 5) 
Missing high redshift clusters in the data used 
(e.g., the EMSS). 6) Evolution of the L-T rela- 
tion. 

Although the quoted uncertainties have so far 
been of minor importance with respect to the 
paucity of observational data, a breakthrough is 
needed in the quality of the theoretical framework 
if high-redshift clusters are to contribute to in the 
high-precision-era of observational cosmology. 

Moreover, the proper comparison of the two- 
point correlation function with X-ray data, re- 
quires better theoretical predictions which take 
account of the selection function of X-ray clusters 
(Kitayama, Sasaki & Suto 1998), the huuinosity- 
and time dependent bias (Mo & White 1996; 
Jing 1998; Moscardini et al. 1998), the light-cone 
effect (Matarrcsc et al. 1997; Matsubara, Suto 
& Szapudi 1997; Nakamura, Matsubara & Suto 
1998; Yamamoto & Suto 1999) and the redshift- 
spacc distortion (Hamilton 1998; Matsubara & 
Suto 1996; Suto et al. 2000; Nishioka & Yamamoto 
1999; Yamamoto, Nishioka & Suto 1999; Magira, 
Jing & Suto 2000). 

The above discussion and recent developments 
in terms of both theory (improved relations for the 
mass function, M-T relation, and bias) and obser- 
vation (X-ray data) suggest that it would be useful 
to recalculate the two-point correlation function 
and to revisit the constraints on cosmological pa- 
rameters obtained until now. 

Likely in Del Popolo (2003), in the present pa- 
per we are principally interested in studying the 
effects of these changes on the values of the cos- 
mological parameters and in comparing them with 
previous estimates, and then in the specific values 
obtained. For this reason, we made a comparison 
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of the theoretical results with observations using 
several samples such as RASSl and XBAC, even 
if it is known that the REFLEX is more precise 
(small errorbars). The paper is organized as fol- 
lows: in Sect. 2, we introduce the model used. Sect. 
3 is devoted to the results and Sect. 4 to discussion 
and conclusions. 



correlation function of clusters in redshift space is 
computed as 

4S(i?,z) = _y^ dkeP!,{k,z)joikR), (3) 
where jo{kR) is the spherical Bessel function. 



2. Theoretical model 

2.1. Redshift-space distortion 

In order to obtain a theoretical model for the spa- 
tial two-point correlation function of galaxies in 
different cosmologies, we follow and improve the 
paper of Suto et al. (2000) (hereafter S2000). Their 
model takes proper account of nonlinear gravi- 
tational evolution of mass fluctuations, redshift- 
space distortion due to the linear peculiar veloc- 
ity field and to finger-of-god effect, cluster abun- 
dance and bias evolution on the basis of the Press 
- Schechter theory, and the light-cone effect. 

As previously reported, one of the effects to 
take into account is the two-point correlation func- 
tion distortions due to the peculiar velocity field. 
We take into account this redshift-space distortion 
following Cole et al. (1994), Magira et al.(2000) 
and Yamamoto et al. (1999). Assuming that the 
bias of the cluster density field relative to the mass 
density field is linear and scale- independent, the 
power spectrum in redshift space is well approxi- 
mated by: 



pS(fc,/i,z)=P«,,3(fc,z) [feel(z)] 



1 + (fc^a^,)2/2 



where P^^^^{k,z) is the mass power spectrum in 
real space, fj, the direction cosine of the wavenum- 
ber vector and the line-of-sight of the fiducial ob- 
server, and (3 is linear redshift-space distortion 
(Kaiser 1987), defined by 



1 d\nDi{z) 
bci{z) dhia{z) ' 



(2) 



2.2. The evolution of the mass auto-correlation 
function 

To predict the clustering properties of X-ray clus- 
ters, we need a description of the matter covari- 
ance function and its redshift evolution. To this 
aim we used the method of Smith et al. (2003), 
which is an improvement of the method of Peacock 
& Dodds (1994), and Peacock & Dodds (1996) 
^ for evolving ^(r, z) into the fully non-linear 
regime. The authors adopted a new approach to 
fitting power spectra, based upon a fusion of the 
halo model and HKLM (Hamilton, Kumar, Lu, 
Matthews) scaling. This approach has been em- 
pirically shown to allow an accurate description of 
a very wide range of power spectrum data. Their 
formula reproduced the scale- free power spectrum 
data and also the CDM results of Jing (1998) with 
an rms error better than 7% (see Smith et al. 
(2003) for more details and their Appendix for the 
fitting formula). 

The linear power spectrum used in this pa- 
per, Pl oc k^^T'^{k), uses the Bardeen et al. 
(1986) transfer function T{k) (Bardeen et al. 
(1986)(Eq. (G3)), and the shape parameter F is 
2 given by: 
(1) 



F 



17om/iexp(-f2ob - y/h/0.5 ilob/i^Om) , (4) 



where bci{z) is the redshift-dependent bias factor 
and Di{z) is the linear growth factor normalized 
to be unity at the present time. 

The denominator in Eq. takes account of 
the nonlinear redshift-space distortion (finger- of- 
God) assuming that the pair-wise velocity distri- 
bution in real space is exponential with the veloc- 
ity dispersion of av{z). 

As in S2000, to calculate we use the fitting 
formula of Mo, Jing & Borner (1997). Averaging 
Eq. (Q) over the angle with respect to the line-of- 
sight of the observer one obtains Pfi{k, z) similarly 
to S2000 (Eq. 4-7). The corresponding two-point 



(Sugiyama 1995), where ilom is the baryonic con- 
tribution to the density parameter. In the part 
of the paper dealing with the direct comparison 
with XBACs and RASSl data, we consider an 
open CDM model (OCDM), with matter density 
parameter Slom — 0.3 and erg = 0.87, and a low- 
density flat CDM model (ACDM), with f^om = 
0.3, and (Tg = 0.93 (see e.g. Liddle et al. 1996a, b 
and references therein). 

In the part of the paper dealing with the 
maximum-likelihood analysis the value of F is al- 
lowed to vary in the range 0.05-0.5, while fiom 
ranges from 0.1 to 1 in the framework of both open 
and flat models. The normalizations of the primor- 
dial power-spectrum, parameterized by erg (the 
r.m.s. fluctuation amplitude in a sphere of 8ft,~^ 
Mpc) is allowed to vary in the range 0.2 < erg < 2. 
In the maximum-likelihood analysis the cosmolog- 
ical models considered, are defined by four param- 
eters: r^om, ^^oA, r and erg. 



Based on Hamilton et al. (1991) original ansatz 
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Fig. 1. The bias factor 
The solid line represents the 
diction of Mo & White (1996) 
prediction for b obtained fro: 
paper and the dashed line ithi 
prediction of Sheth & Tormeri 



as a i function of . v^ 



^^|herical collapse pre- 
the dotted line the 
the model of this 
^ ellipsoidal collapse 
(1999). 






Fig. 2 

Fig. 2. Comparison of various mass functions. The 
dotted line represents the Sheth & Tormen (2002) 
prediction, the solid line that of Jenkins et al. 
(2001) and the dashed line that of Del Popolo 
(2002b). 



2.3. Bias evolution 

In order to predict the clustering properties 
of clusters as a function of redshift, we need 
to know how bias evolves. S2000 adopted for 
the 'monochromatic' bias b{M, z) the expression 
which holds for virialized dark matter haloes (see 
their Eq. 17) (e.g. Mo & White 1996; Catelan et al. 

1998) , and to get the effective bias factor the Mo 
& White (1996) equation was combined with the 
Press-Schechter relation to translate the quoted 
bias factor into a function of X-ray flux limit (sec 
their Eq. 18). 

Several papers in literature has shown that the 
Mo & White (1996) bias formula does not cor- 
rectly reproduce the correlation of low mass haloes 
in numerical simulations. 

Several alternative fits have been proposed 
(Del Popolo & Gambera 1998; Jing 1998; Porciani, 
Catelan & Lacey 1999; Sheth & Tormen 1999; Jing 
1999; Del Popolo 2001). The bias model of Sheth 
& Tormen (1999) has been shown to produce an 
accurate fit of the distribution of the halo popu- 
lations in the GIF simulations (Kauffmann et al. 

1999) . 



(b) 

(a) M-T relation predicted by the modified 
continous cluster forr^ation model. The solid line 
is the prediction of the: quo ted model, described in 
Del Popolo '(2662a)), for jIa = 0, f^o = 1, shifted 
downwards, similarly as in Afshordi & Cen (2002), 
to fit the Finoguenov et al. (2001) observational 
data at the massive end. The dotted line repre- 
sents the prediction of the same model for f^A = 
0.7 and f^o — 0.3. (b) Temperature evolution pre- 
dicted by the modified continous cluster formation 
model in the case fio — 0.3. The dotted line repre- 
sents the "classical" prediction, Tx « [l + z). The 
short-dashed line represents the late-formation 
approximation as expressed by Eq. (8) in Voit 

(2000), namely Tx cx [-^] ^'\^ + z)- The 
long-dashed line and the dot-dashed one plot Eq. 
(40) in Del Popolo (2002a) for the spherical col- 
lapse model n=— 2,n = — 1, respectively. The 
solid line plots the same equation for n — —\ and 
taking account of angular momentum acquisition 
by the protostructure. 

In this paper we adopt the bias model de- 
scribed in Del Popolo & Gambera (1998), and 
Del Popolo (2001), because it produces a mass 
function that is in better agreement with the pre- 
dictions of Jenkins et al. (2001) (see below). The 
quoted biased model is based on the threshold: 



1 



Tte.L'^ ■ dr 



(5) 



where 5co — 1-68 is the critical threshold for a 
spherical model, n is the initial radius, rta is the 
turn-around radius, and L is the angular momen- 
tum. ^ Eq. 0, is obtained taking account of the 

The angular momentum appearing in Eq. (|KJ is 
the total angular momentum acquired by the proto- 
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effect of asphericities. In Fig. 1, we plot the bias 
parameter, b, as a function of the peak height v, ^ 
which is proportional to the halo mass. The solid 
line shows the spherical collapse prediction of Mo 
& White (1996), the dotted line the prediction for 
h obtained from our model, and the dashed line the 
ellipsoidal collapse prediction of Sheth & Tormen 
(1999). As shown in the figure, taking account 
of the effects of asphericity produces a change in 
the dependence of 6 on i/ in good agreement with 
Sheth & Tormen (1999). Prom Fig. 1 it is evi- 
dent that at the low mass end the bias relation 
has an upturn, meaning that less massive haloes 
are more strongly clustered than in the prediction 
of the spherical collapse model of Mo & White 
(1996) and in agreement with N-body simulations 
(Jing 1998; Sheth & Lemson 1999; Kauffmann et 
ah 1999). 

The above bias factor can be translated into a 
function of X-ray fltix limit according to: 



Assuming that the barrier is proportional to 
the threshold for the collapse, as in ST, the barrier 
can be expressed in the form: 



B{M) 



ben{z, > Slim) = 



Mu„(Sn„) 



dM b{z,M) n{z,M) 



/ dM n{z, M) 

JMum 



-(6) 



where n{z,M) is the number of objects actually 
present in the catalogue with redshift in the range 
z, z+dz and M in the range M, M+dM. One can 
estimate n{z, M) from the Press-Schechter (1974) 
(PS) formula; however, several studies have shown 
some discrepancies between PS and simulations. 

In order to obtain a better estimate of n(M, z), 
wc can use the excursion set approach, which al- 
lows one to calculate good approximations to sev- 
eral important quantities, such as the "uncondi- 
tional" and "conditional" mass functions. Sheth 
& Tormen (2002) (hereafter ST) provided formu- 
las to calculate these last quantities starting from 
the shape of the barrier. They also showed that 
the "unconditional" mass function, which is the 
one we need now, is in good agreement with re- 
sults from numerical simulations. Using the bar- 
rier shape obtained in Del Popolo & Gambera 
(1998), obtained from the parameterization of the 
nonlinear collapse discussed in that paper, to- 
gether with the results of ST we can calculate the 
"unconditional" mass function. 



structure during evolution. In order to calculate L, We 
will use the same model as described in Del Popolo & 
Gambera (1998, 1999) (more hints on the model and 
some of the model limits can be found in Del Popolo, 
Ercan & Gambera (2001)). 

1/ = ^ ct(m) ) ^® ratio between the critical over- 
density required for collapse in the spherical model, 
Sc{z), to the r.m.s. density fluctuation a{M), on the 
scale r of the initial size of the object M 



1 + 



rta,L^ ■ dr 



1 



(7) 



(Del Popolo 2001), where ^co = 1-68 is the critical 

threshold for a spherical model, r; is the initial 
radius, rta is the turn-around radius, L the angular 
momentum, a\ = 0.585 and /3i = 0.46. 

As described in Del Popolo (2002b), the mass 
function can be approximated by: 



n{M, z) 



p d log V 
M2 dlogM 



yf{u) ~ 1.21 




M2 dlogM 



where a = 0.707, and p is the background density 
^. In Fig. 2 we plot a comparison of the various 
mass functions: the dotted line represents Sheth 
& Tormen (2002) prediction, the solid line that 
of Jenkins et al. (2001) and the dashed line that 
of Del Popolo (2002b). As Fig. 2 shows, the mass 
function obtained in this paper is in very good 
agreement with that of Jenkins et al. (2001) in the 
regime probed by the simulations. Notice that a 
large part of the constraints obtained in the past 
has been obtained using the PS mass function. 
Only in some more recent papers has the mass 
function been calculated by means of ST model 
(e.g. Borgani et al. 2001) or that of Jenkins et al. 
(2001) (Hamana et al 2001; MMM). Moreover the 
M-T relation chosen is the usual one obtained sim- 
ply from the virial theorem (see next subsection) . 
In other words, this paper introduces noteworthy 
improvements on the previous calculations in lit- 
erature. 

Before concluding this subsection, following 
Hamana et al. (2001), we further attempt to in- 
corporate the scale-dependence on the basis of the 
results of Taruya fc Suto (2000), in which the 
scale-dependence arises as a natural consequence 
of the formation epoch distribution of halos. 
Yoshikawa et al. (2001), had shown that the scale- 
dependence of the Taruya & Suto (2000) model 
agrees with their numerical simulations. Therefore 
we construct an empirical halo bias model of the 
two-point statistics which reproduces the scale- 
dependence of the Taruya & Suto (2000) bias with 
the amplitude fixed by the mass-dependent bias 
obtained in Del Popolo & Gambera (1998) and Del 
Popolo (2001). Following Hamana et al. (2001), 
the scale dependent bias shall be described by the 
following simple fitting formula: 

6h(M, R, z) = b{M, z) [1.0 + b{M, z)aR{R, z)f^^ , (9) 

2 



Note that in this formula u ■ 



( icCfiV 

\a(M)) 
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for R > 2Rvi,:{M,z), and otherwise 0, where 
Rvir{M,z) is the virial radius of the halo of 
mass M at z and <tr[R^z) is the mass variance 
smoothed over the top-hat radius R. The bias fac- 
tor b contained in Eq. Q through the spectrum 
P^^ should then be substituted by Eq. Q . 

Although the modeling is not completely self- 
consistent in the sense that the scale-dependence 
of the halo biasing factor is neglected in describ- 
ing the redshift distortion, the above prescription 
provides a good approximation since the scale- 
dependence in the biasing is of secondary impor- 
tance in the redshift distortion effect of halos (see 
Hamana et al. 2001). 

2.4. Limiting flux and halo mass 

In order to predict the abundance and clustering 
of X-ray clusters in a given sample, it is neces- 
sary to relate the X-ray fluxes to the correspond- 
ing halo mass at each redshift. As a first step, we 
relate the total mass of the dark halo of a cluster 
to the temperature of the gas. The M-T relation 
that we shall use is that calculated in Del Popolo 
(2002a). The M-T relation is calculated using the 
merging-halo formalism of Lacey & Cole (1993), 
which takes account of the fact that massive clus- 
ters accrete matter quasi-continuously, and the 
present paper is an improvement of a model pro- 
posed by Voit (2000) (hereafter V2000), again to 
take account of angular momentum acquisition by 
protostructures and of an external pressure term 
in the virial theorem. The M-T relation obtained 
in Del Popolo (2002a), is given by: 



is a break at the low mass end (T ~ 3 — 4keV) of 
the M-T relation is present. The behavior of the 
M-T relation is as usual, M oc T^/^, at the high 
mass end, and M oc T'^ , with a value of 7 > 3/2 
dependening on the chosen cosmology. Larger val- 
ues of 7 are related to open cosmologies, while 
ACDM cosmologies give a slope intermediate be- 
tween the flat case and the open case. 

In Fig. (3a), we plot the M-T relation predicted 
by the modified continous formation model. The 
solid line is the prediction of the quoted model, 
described in Del Popolo (2002a), for Vt/v = 0, 
r^o = Ij shifted downwards, as in Afshordi & 
Cen (2002), to fit the FRB observational data 
at the massive end. The dotted line represents 
the prediction of the same model for JIa = 0.7 
and VIq = 0.3. In Fig. (3b), we plot the temper- 
ature evolution predicted by the modified model 
for r^o — 0.3. The dotted line represents the "clas- 
sical" prediction, Tx oc (1 -I- z). The short-dashed 
line represents the late-formation approximation 
as expressed by Eq. (8) in Voit (2000), namely 



kT ~ SkeV 



Aft 



3 I K{m,x) 



MS/3 



10i5/i-iM( 



o 



where ta 



„ 3 , m 

•noA)2 



5/{n + 3) ( 

Ho(l-no™-aoA)2 

where n is the spectral index), Mq is defined in 
Del Popolo (2002a), and: 

K{m,x) = Fx (m — 1) LerchPhi{x,l,3m/5 + 1) — 
F{m-l)LerchPhi{x,l,3m/5), (11) 

where F is defined in Del Popolo (2002a) (Eq. 35) 
and the LerchPhi function is defined as follows: 



LerchPhi{z, a, v) — 



n—0 ^ ' 



(12) 



and where KQ{m,x) indicates that K{m,x) must 
be calculated assuming t — t^. 

Eq. (|10ll takes account of the fact that massive 
clusters accrete matter quasi-continuously, and 
also of tidal interaction between clusters. The ob- 
tained M-T relation is no longer self-similar, there 



Tx oc A. 



1/3 



1/3 



(l + z). The long-dashed line 

and the dot-dashed one plot Eq. (40) in Del Popolo 
(2002a) for the spherical collapse model n = —2, 
n — ~1, respectively. The solid line plots the same 
equation for n — —1 and taking account of angular 
momentum acquisition by the protostructure. 

The next step (see S2000) is to transform the 
temperature to the luminosity of clusters using 
the observed luminosity-temperature relation. In 
S2000, they assumed: 







Ko{m,x) 









Lhol ~ L. 



44 



( n 



V 6keV / 



(l+z)« lO^^h- 



^erg sec ^(13) 



^^itWL^^ = 2.9, a = 3.4 and C = on the basis of 
the quoted observational indications (e.g., David 
et al. 1993; Ebeling et al. 1996; Ponman et al. 
1996; Mushotzky & Scharf 1997). 

Several independent analyses of nearby clus- 
ters with Tx > 1 keV consistently show that 
L44 ~ 3 and a ~ 2.5 — 3.5 (e.g., White, Jones 
& Forman 1997; Wu, Xue & Fang 1999, and ref- 
erences therein). For cooler groups, < 1 keV, the 
Lho\ — T relation steepens. Mushotzky & Scharf 
(1997) found that data out to z ~ 0.4 are consis- 
tent with no redshift evolution in the Lboi — T re- 
lation out to z ~ 0.4. In Moscardini et al. (2000a, 
2000b) the authors translated the cluster bolomet- 
ric luminosity into a temperature, adopting the 
empirical relation 



(14) 



where the temperature is expressed in keV and 
Lboi is in units of 10^'' erg/s and A = 4.2 and 
B = 1/3, which are a good representation of the 
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data with T > 1 keV (e.g. Markevitch 1998 and 
references therein). 

From what was previously said, it is clear that 
the L-T relation is a source of uncertainties. As 
in Del Popolo (2003), in the present paper we 
are principally interested in studying the effects 
of the improvements on the M-T relation, mass 
function and bias model on cosmological param- 
eters and to compare them with previous esti- 
mates. For this reason, in the following, we shall 
follow the philosophy of Moscardini et al. (2000a), 
Borgani, Plionis, & Kolokotronis (1999), namely 
we shall adopt a 'default' value for a and L44, 
(L44 = 2.9, and a — 3.4 as the reference values), 
and we calculate the correlation function and the 
constraints on cosmological parameters. We shall 
compare these with the results of previous pa- 
pers (e.g., Moscardini et al. 2000a, 2000b; S2000). 
Finally, we shall calculate the effects of the vari- 
ation of OL (in the range 2.5 < a < 3.5) on the 
resulting model constraints. Notice that all plots 
shown in the next sections are based on these 'de- 
fault' values L44 = 2.9, and a — 3.4. 

After fixing the L-T relation, iboiCTgas) 
is transformed in the band-limited luminosity 
-f'band[Tgas,£^i,£^2] as showu by S2000 (Sect. 2.2). 

To obtain il/ii„i, necessary to calculate 6cff in 
Eq. @, we use the method of S2000 (see their 
Sect. 2.4) 

2.5. The light-cone effect 

The final step is to calculate the two-point cor- 
relation function on the light cone (Yamamoto & 
Suto 1999): 

^X-cl(^; > 'S'lim) — ^Tz— ^ 

dz — - n'i.{z) 
dz 

where R is the comoving separation of a pair of 
clusters, Zmax and Zmin denote the redshift range of 
the survey, and £,^\{R, z; > S'lim) is the correspond- 
ing two-point correlation function on a constant- 
time hypersurface at z in redshift space (Eq.|(2Jl). 
The comoving number density of clusters in the 
flux-limited survey, uqIz; > S'lim), is computed by 
integrating the mass function Eq. © as: 

/•OO 

710(2;; > Slim) = / n{M,z)dM. 

Finally the comoving volume element per unit 
solid angle is 

dVc _ dljz) 

m ~ H{z) ■ 

where: 



2.6. Maximum-likelihood analysis 

In order to obtain constraints for cosmological 
models, we use a maximum-likelihood analysis. 
One possibility to accomplish the quoted analysis 
is as shown by Marshall et al. (1983), Del Popolo 
(2003), or Borgani, Plionis & Kolokotronis (1998). 

In the present paper. We used the same model 
as MMM. The likelihood is C oc exp(-x^/2), 
where 



X 



[ro(z) -ro(i;f}om,r^oA,r, as)]" 



(16) 



The sum runs over the observational dataset de- 
scribed in Sect. 2 of MMM, i.e. A^data = 3 and 
A'data = 4 for the optical and X-ray bands, respec- 
tively (in the present paper we shall use only X-ray 
data). The quantities ro(i) and (Jro{i) represent 
the values of the correlation length and its la er- 
rorbar for each catalogue, as reported in Table 1 of 
MMM; ro(i; ilom, ^^OA, r, (Tg) is the corresponding 
theoretical prediction obtained with a given choice 
of cosmological parameters. The best-fit cosmo- 
logical parameters are obtained by maximizing £, 
i.e. by minimizing . The 95.4 and 99.73 per cent 
confidence levels for the parameters are computed 
by finding the region corresponding to an increase 
A^2 with respect to the minimum value of x^. 
Other details of the maximum-likelihood analysis 
are given in the next section. 

3. Results 

We begin the analysis by comparing the theoreti- 
cal predictions for the two-point correlation func- 
tion to the observational clustering properties of 
!iiA:^Sl and afterwards to those of XBACs. 
Fig. 4 shows the comparison with RASSl. In 



the plot, the dashed line represents the predic- 
tion of Moscardini et al. (2000b) for their ACDM 
model, and the solid line represents the predic- 
tion of the ACDM, with f^om = 0.3, erg = 0.93, 
and r = 0.21 (see e.g. Liddle et al. 1996a,b 
and references therein), for the 'default' value of 
a = 3.4. The plot for the OCDM model was not 
reported since it gives very similar predictions to 
the ACDM model. The ACDM model is in good 
agreement with the data and the predictions are 
always inside the 1 tr errorbars, tq ~ 28.2h~^ Mpc. 
A comparison with the Moscardini et al. (2000b) 
results, after the latter has been corrected taking 
account of the description of clustering in the past- 
light cone, (which gives ~ 22. 4h^^ Mpc), shows 
that in our model the correlations are higher, by 
more than 20%. Varying the value of a in the 
range 2.5-3.5, not plotted in Fig. 4, the result- 
ing spa tial correlation function does not change 



H{z) - i?oVf^Om(l + zY + (1 - f^om - fioA)(l + 2)it^ fMa(|)(i,^ut Only its amplitude: the smaller a is, 
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Fig. 4. Comparison of the cluster spaee eorre- 
lations in the RASSl sample with the theor^iO; 
model of the present paper. In the plot, the dashed 
line represents the prediction of Moscardini et 
(2000b) for their ACDM model, and the solid 
that of the ACDM (l^om = 0.3, as = 
lated using the model of this paper. 



Fig. 5 




r(Mpc/h) 



5^ 
-14 



-13 



-12 



-11 



Fig. 5. Comparison of the cluster space correla- 
tions in the XBACs sample with the theoretical 
model of this paper. The observational estimates 
are shown by two regions: the first (enclosed in the 
solid lines connected by the vertical solid lines), 
refers to the (Icr) estimates obtained by Abadi, 
Lambas & Muriel (1998), while the second region 
(enclosed in the dashed lines connected by the hor- 
izontal dashed lines), shows the (2(7) estimates by 
Borgani, Plionis & Kolokotronis (1999). The solid 
curve represents the Clom = 0.3, erg = 0.93 ACDM 
model calculated using the model of this paper, 
while the dashed line represents the prediction of 
Moscardini et al. (2000a) for their ACDM model. 



the smaller ^(r) is. However, the changes are quite 

small, as shown in Fig. 11 of Moscardini et al. 
(2000a). The induced change of ro, can be written 
as ro 28.2 ± 5.2h-^ Mpc. 

In Fig. 5 the model is compared with the 
XBACs catalogue. Abadi, Lambas & Muriel 

(1998) found that the XBACs spatial correla- 
tion function can be fitted by the usual power- 
law relation ^(r) = {r/ro)~^ with 7 = 1.92 and 
ro = 21. 1+2 3 ^P"^ errors correspond 
to 1 cr), while Borgani, Plionis & Kolokotronis 

(1999) , who adopted an analytical approximation 
to the bootstrap errors for the variance of ^(r), 
found 7 = 1.98t^;i and a slightly larger value of 



log S,^ (erg s ' cm ^) 



Fig. 6. The behavior of the correlation length ro 
as a function of the limiting X-ray flux Sum ■ The 
solid line represents the ACDM model according 
to the model of this paper, and the dashed line 
that of Moscardini et al. (2000a) 



2-C7 uncertainties). In Fig. 5, these observational 
estimates are compared with the theoretical pre- 
dictions of this paper. The observational estimates 
are shown by two regions: the first one (enclosed 
in the solid lines connected by the vertical solid 
lines), refers to the (1 a) estimates obtained by 
Abadi, Lambas & Muriel (1998), while the second 
region (enclosed in the dashed lines connected by 
the horizontal dashed lines), shows the (2 a) esti- 
mates by Borgani, Plionis & Kolokotronis (1999). 
The solid curve represents the ilom = 0.3, erg = 
0.93 ACDM model, while the dashed line repre- 
sents the prediction of Moscardini et al. (2000a) 
for their ACDM model 

The results are in qualitative agreement with 
MMM, but the value for ro (ro ~ 28.2h-i Mpc) 
obtained here is larger than the value in the 
quoted paper (ro ~ 22h~'^ Mpc for ACDM), while 
it is in better agreement with that of Borgani, 



ro = 26.0t4;7 ^^"^ (the errors in this case are 28.2h"^ Mpc). 



® For the ACDM model (solid line), the results are 
in good agreement with the observational data (ro ^ 
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Fig. 8. 




Confidence contours of F and ilom- 
Dashed and solid lines represent 95.4 and 99.73 
per ccn; confidence levels. The left panel refers to 
flat cosmological models with varying cosmolog- 
ical cor st ant the right one to open models with 
vanishiiig Oqa. 



0.2 




0.2 0.4 



0.4 0.6 

n 



(b) 



Fig. 7. The variation of A 



fitting value of the matter 
for flat CDM models (with 
constant Qq\ = 1 — fiom) with 
r = 0.2 and normalization re 
ter abundance. The solid line 
of MMM (obtained using the X 
ter dataset), and the dashed 
per. Horizontal lines correspond^^ng 
99.73 per cent confidence levels 




around the best- 
density parameter fiflij 
varying cosmologieal 
shape parameter 
ucing /the clus- 
rebresents /the result/ 
fay complete clus-' 
that iof this pa 
to the9^.4a;nd 
are also shown. 



pjiodr 



lir[e 



Plionis & Kolokotronis (1999), who found a value 
of ro = 26l4;7h^^ Mpc. Similar coi^gideratidijs to 
those relative to Fig. 4 are valid if we vary the 
value of a. 

In order to study the possible dependence of 
the clustering properties of the X-ray clusters on 
the observational characteristics defining the sur- 
vey, we plot the values of the correlation length 
ro in the catalogues where we vary the limiting 
X-ray flux S'lim- The result is shown in Fig. 6: the 
dashed line is the result obtained by MMM, and 
the solid line that using the model of this paper, in 
the case of a ACDM model. The other cosmolog- 
ical models, not plotted, have a similar behavior 
but smaller amplitude (see MMM). In Fig. 6, ro 
increases with Luni- The increase is more rapid 
for ACDM models (see MMM). A similar anal- 



(b) 

Fig. 9. Confidence contours of F and as. Dotted 
and solid lines represent 95.4 and 99.73 per cent 
confidence levels. Panel (a) refers to a CDM model 
with riom = 0.3, QoA = 0. Panel (b) refers to a 
CDM model with Qq^ = 0.3, Ooa = 0.7. 



ysis has been made by S2000 (see their Fig. 8) 
and MMM. Even if the S2000 results cannot be 
directly compared with those of MMM, (differ- 
ences in cosmological parameters, formalism for 
the past-light cone effect, etc.) or even with those 
of this paper (differences in cosmological models 



10 



A. Del Popolo et al.: Cosmological constraints from galaxy clusters. 



parameters adopted), it is clear that there is a 
qualitative agreement but MMM tends to predict 
smaller correlation lengths (by approximately 20- 
30 %) with respect to S2000 and the latter pre- 
dicts smaller correlation lengths than this paper. 
Part of the difference between MMM and S2000 
comes from the different values of the exponent 
of the temperature used by the quoted authors 
in the temperature-luminosity relation: in MMM 
B = 1/3 (see Sect. 2.4)), and that used in S2000 
{a = 1/3.4 (see their Eq. 11)). Moreover, in their 
approach S2000 include a method to take account 
of redshift-space distortion effects not completely 
considered by MMM which tends to increase the 
correlation estimates. 

In the present paper, the correlation lengths 
are larger in comparison with both the previous 
two quoted papers predictions. These differences, 
expecially when we compare the results with those 
of S2000, are due essentially to the different M- 
T relation, mass function and bias model used 
in this paper. As reported in Sect. (2), the fun- 
damental goal of this paper is to study the ef- 
fects of the improvements quoted above on the 
values of the, cosmological parameters, and then 
to constrain the cosmological parameters by the 
clustering properties of clusters of galaxy using a 
maximum-likelihood analysis. This analysis was 
started by considering as a free parameter flm 
only, fixing T = 0.2, which is in the range sug- 
gested by various other works (see e.g. Peacock 
& Dodds 1996), and ag to reproduce the cluster 
abundance. 

For the normalization we adopt the fitting for- 
mula by Pierpaoli et al. (2003), with the M-T nor- 
malization parameter T* = 1.75 (see Pierpaoli et 
al. 2003). 

In Fig. 7, we plot the results of the maxi- 
mum likelihood analysis obtained by using X-ray 
band data (RASSl, BCS, XBACs, REFLEX) for 
flat models with varying cosmological constant 
(OoA = 1 - i^om), with F = 0.2. The solid line 
represents the results obtained using the complete 
X-ray data set by MMM, and the dotted line gives 
the results for the model in the present paper. 

The constraints on fim are: 0.25 < fJom < 0.45 
and 0.23 < l^om < 0.52 at 95.4 and 99.73 per cent 
levels, respectively. The result shows larger values 
for riom on the order of 20-30 %, when compared 
with MMM results (0.2 < f^om < 0.35 and 0.2 < 
Oom < 0.45 at the 95.4 and 99.73 per cent levels, 
respectively) . 

In Fig. 8 we show the results of the maxi- 
mum likelihood analysis fixing the model normal- 
ization to rejiToduce the cluster abundance, and 

They used a zero-order model similar to that of 
Kaiser 



leaving two free parameters: F, and f^mo- The solid 
lines represents the 99.73 confidence levels, and 
the dashed lines the 95.4 confidence levels. The 
figure shows that the allowed regions, at least for 
fiom > 0.5, depend strongly on F. In the case of 
the fiat model (Fig. (8a)), the constraints for F are 
0.1 < F < 0.14, and in the case of the open model 
(Fig. (8b)) 0.09 < r < 0.13. These values are 
smaller than those of MMM, by about 20 - 30%. 
(Note that in MMM paper, the reported values 
are approximated. A comparison of our and their 
confidence contours shows the difference). 

In Fig. 9 we show the constraints in the T — 
plane (after keeping the values of J^a fixed). We 
consider an open model with ilom = 0.3, Fig. (9a), 
and a flat model. Fig. (9b), again with ^^om = 0.3. 

8 

For the open model with Oom = 0.3 the 2a 
region has F in the range 0.11-0.2 and as between 
0.7 and 1.55. The values of as obtained here are 
larger by ~ 20% than those of MMM. 

For a fiat model with fiom = 0.3 the 2a region 
has 0.13 < F < 0.2 and 0.8 < as < 1.6. The values 
of as obtained are larger than those of MMM by 
-20% 

The effect of varying a on the maximum- 
likelihood analysis has been also studied. Prom the 
study, it turns out that F is quite insensitive to a 
change of a, while erg is only weakly dependent on 
it: the minimum in the maximum-likelihood anal- 
ysis decreases (increases) by ~ 10% for a = 2.5, 
(for a = 3.5) (see also Fig. 8 of Borgani, Plionis & 
Kolokotronis 1999). For Oom the changes are sim- 
ilar to those of F but larger (of the order of 20 %) 
(see also Fig. 4 of Borgani et al. 2001, and Fig. 10 
of Eke et al. 1998). 

Notice that in the analysis of this paper, we 
used only the data from the X-ray catalogues and 
not a combination of optical and X-ray catalogues. 
This is because, as shown by MMM, the combina- 
tion of optical and X-ray catalogues gives results 
that are almost indistinguishable from those ob- 
tained by the X-ray analysis only. As in MMM, 
the constraints from the X-ray datasets are in gen- 
eral tighter than those obtained from the optical 
data. 

Although the quoted uncertainties have so 
far been of minor importance because of the 
paucity of observational data, a breakthrough is 
needed in the quality of the theoretical model- 
ing if high-redshift clusters are to take part in the 
high-precision-era of observational cosmology. As 
shown, using models like PS for the mass func- 
tion instead of models that are in agreement with 

* Note that the value of Qom chosen is not too 
different from that obtained by WMAP: ^omh^ = 
0.14 ± 0.02 and h = 0.72 ± 0.05 (Spergel et al. 2003) 
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Jenkins et al. (2001), introduces errors of the or- 
der of 20-30 % in the values of the parameters 
constrained (see Del Popolo 2003) and other er- 
rors are introduced if one uses a simphfied ver- 
sion of the M-T relation. Moreover a proper in- 
terpretation of such redshift surveys in terms of 
the clustering evolution requires an understand- 
ing of many cosmological effects which can be ne- 
glected for << 1 and thus have not been consid- 
ered seriously so far. For example, neglecting the 
light-cone effect leads to underestimates of ~ 20% 
in samples like RASSl, and of up to 25 % for a 
deeper survey such ABRIXAS (Moscardini et al. 
2000a) . Neglecting red-shift space distortions pro- 
duces underestimates of ro of ~ 10%: as reported 
in the present paper, a comparison between the 
behavior of the correlation length as a function 
of the limiting X-ray flux, in papers like S2000 
and MMM shows a difference of ~ 30% because 
in MMM the redshift-space distortion effects were 
not taken into account, and because of a difference 
in the exponent in the temperature-luminosity re- 
lation. Taking into account the asphericity in grav- 
itational collapse that leads to a different relation 
for bias (Del Popolo & Gambera 1998; Sheth & 
Tormen 1999), to a modified version of the mass 
function (Del Popolo & Gambera 1998; Sheth & 
Tormen 2002; Del Popolo 2002b), and a different 
M-T relation (Del Popolo 2002a), leads to higher 
values of Oom, erg, fo by at least 20 %. In the near 
future theory and observations should converge to- 
wards a more precise constraining of cosmological 
parameters. 

4. Discussion and Conclusions 

In this paper, we have recalculated the two-point 
correlation function of clusters of galaxies for 
OCDM and ACDM cosmological models, improv- 
ing the model of S2000, in the light of recent 
theoretical developments, by using the theoreti- 
cal mass function derived in Del Popolo (2002a), 
the M-T relation derived in Del Popolo (2002b) 
and the bias model of Del Popolo (2001). As in 
Suto's paper, the model properly takes accoimt of 
nonlinear gravitational evolution of mass fluctu- 
ations, redshift-space distortion due to the linear 
peculiar velocity field and to fingcr-of-god, cluster 
abundance and bias evolution, and the light-cone 
effect. This theoretical model has before been com- 
pared with the observed spatial correlation func- 
tion for clusters in RASSl, and in XBACs sam- 
ples. The comparison shows that only the pre- 
dictions of models with = 0.3 arc in good 
agreement with data. The results are in qualita- 
tive agreement with MMM, but the values for ro 
here obtained (ro ~ 28.2h-i Mpc for the ACDM 
model) are larger than the values of the quoted pa- 



per (MMM), (ro ~ 22h-i Mpc for ACDM), while 
they are in better agreement with that of Borgani, 
Plionis & Kolokotronis (1999), who found a value 
of ro = 26t4.7h~^ Mpc. In order to study the pos- 
sible dependence of the clustering properties of the 
X-ray clusters on the observational characteristics 
defining the survey, we plot the values of the cor- 
relation length ro in the catalogues where we vary 
the limiting X-ray flux Sum- ^ All the cosmolog- 
ical models displays an increase of ro with Liim- 
The increase is more rapid for ACDM models (see 
MMM). Comparing the result with those of a sim- 
ilar analysis by S2000, (see their Fig. 8) and MMM 
it is clear that there is a qualitative agreement but 
MMM tends to predict smaller correlation lengths 
(by approximately 20-30 %) with respect to S2000, 
and the latter predicts smaller correlation lengths 
than the present paper. These differences, expe- 
cially when we compare the results with those of 
S2000, are due essentially to the different M-T re- 
lation and mass function used in this paper. In or- 
der to obtain constraints on cosmological parame- 
ters we performed a maximum-likelihood analysis 
by comparing the theoretical predictions to a set 
of observational data in the X-ray band (RASSl 
Bright Sample, BCS, XBACs, REFLEX), simi- 
larly to MMM. The parameters to be constrained 
are: fim, ^a, the power-spectrum shape parame- 
ter F and the normalization <Js- The constraints 
obtained for the matter density parameter in a 
flat CDM model, are: 0.25 < l^om < 0.45 and 
0.23 < l^om < 0.52 at the 95.4 and 99.73 per 
cent levels, respectively, larger by ~ 20% than 
the MMM predictions. Keeping the model normal- 
ization fixed to reproduce the cluster abundance, 
and leaving two free parameters: F, and flmo, we 
find that for the flat model the constraints for F 
are 0.1 < F < 0.14, while for the open model 
0.09 < F < 0.13. These values are smaller than 
those of MMM, by about 20 - 30%. After flxing 
the values of Oa, we obtain the constraints in the 
F — o"8 plane, showing that if we keep the value of 
l^A fixed the open model with Oom = 0.3 the 2a 
region has F in the range 0.11-0.2 and ag between 
0.7 and 1.55. In the flat model with f^om = 0.3 the 
2a region has 0.13 < F < 0.2 and 0.8 < cts < 1.6 
In all three cases, the values of ag obtained are 
larger than those of MMM by ~ 20%. Varying 
a, it turns out that F is quite insensitive to the 
change of a, while ag is only weakly dependent on 
it: the minimum in the maximum- likelihood anal- 
ysis decreases (increases) by ~ 10% for a = 2.5, 



Notice that this analysis can be related to the 
study of the richness dependence of the cluster corre- 
lation function. In fact, a change in the observational 
limits implies a change in the expected mean inter- 
cluster separation dc- 
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(for a = 3.5). In the case of ilom the changes are 
similar to those of T but larger (of the order of 
20 %). Allowing the shape parameter to vary, we 
find that the clustering properties of clusters are 
almost independent of the matter density parame- 
ter and of the presence of a cosmological constant, 
while they appear to be strongly dependent on the 
shape parameter. The constraints from X-ray data 
are tighter than those coming from optical data. 
In conclusion, the data on clustering properties of 
galaxies can be used to constrain important cos- 
mological parameters like fiom, T and erg. 
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